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Abstract. The primordial density fluctuation inevitably couples to all forms of 
matter via loop corrections and depends on the ambient conditions while inflation 
was ongoing. This gives us an opportunity to observe processes which were in progress 
while the universe was inflating, provided they were sufficiently dramatic to overcome 
suppression by powers of (H/Mp) 2 as 10~ 9 , where H is the Hubble scale during 
inflation and Mp is the Planck mass. As an example, if a primordial magnetic field was 
synthesized during inflation, as suggested by some interpretations of the apparently 
universal 10 -6 gauss field observed on galactic scales, then this could leave traces 
in inflationary observables. In this paper, I compute corrections to the spectrum and 
bispectrum generated by a varying electromagnetic coupling during inflation, assuming 
that the variation in this coupling is mediated by interaction with a collection of light 
scalar fields. If the mass scale associated with this interaction is too far below the 
Planck scale then the stability of perturbation theory can be upset. For the mass-scale 
which is relevant in the standard magnetogenesis scenario, however, the theory is stable 
and the model is apparently consistent with observational constraints. 

Keywords: Inflation, Cosmological perturbation theory, Physics of the early 
universe, Quantum field theory in curved spacetime. 
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1. Introduction 

Our understanding of the very early universe has settled down over the last several 
years to give, on balance, a broadly consistent picture. In this concordance model, an 
adiabatic perturbation is synthesized on superhorizon scales during an early epoch of 
inflation and is converted by gravitational collapse into the observed distribution of 
baryons and cold dark matter in the universe. It is generally accepted that this simple 
timeline is sufficient to explain the vast majority of observational data. 

A zoo of optional components can be added to improve the fit to certain small 
anomalies in the data, or to provide an origin for features of our macroscopic world 
which are presently without a theoretical rationale. Examples of these optional features 
are scalar isocurvature modes (which may source non-gaussianities in the temperature 
anistropy [HE]), vector modes (which may source statistical anisotropy [31 HI [5]), 
or tensor gravitational waves. All of these must be subdominant to the adiabatic 
perturbation. We usually assume that these optional extras can be added or taken 
away without penalty. However, if they mediate sufficiently dramatic processes then we 
must remember that the observable adiabatic fluctuation will couple to all of them 
via loop corrections. In this paper the loop corrections arising from an especially 
interesting optional addition are studied, namely a varying electromagnetic^ coupling 
which is mediated by interaction with a collection of light scalar fields. Among other 
effects, a varying coupling of this sort could be responsible for generating a primordial 
magnetic field. Such magnetic fields are required by observation [7j and although a 
number of possible formation mechanisms are known — see, for example, Ref. [Bj for 
a non- inflationary example — their large-scale coherence and homogeneity hints at an 
inflationary origin parallel with the origin of the primordial density fluctuation. 

What can we hope to learn from the study of such loop effects? One motivation 
is that, despite increasing high-quality experimental data which suggests the simplest 
predictions of inflation are a good fit for observation, it is still unclear whether inflation 
is the right model of microphysics or simply a good parametrization of an approximately 
gaussian, scale-invariant and adiabatic spectrum. One way to approach this question 

f In fact, the coupling in question cannot be the electromagnetic coupling a, the fine structure 
"constant" , because inflation is usually supposed to take place at an energy scale far above the scale 
of electroweak symmetry breaking, and at such energies the electromagnetic field has not yet obtained 
a separate identity. 

It is known that non-Abelian gauge fields are shielded from obtaining a perturbation during inflation 
by their strong self-interactions [6], and therefore the gauge coupling in question must be the U(l) 
Standard Model hypercharge field. A fluctuation in this field will communicate some perturbation 
to the electromagnetic field at electroweak symmetry breaking. However, the detailed identity of the 
Abelian gauge field in question is irrelevant for the question posed in the present paper. For simplicity, 
I will refer to this field as the "electromagnetic" field and its coupling as the "electromagnetic coupling" 
throughout, with the understanding that for applications to magnetogenesis the hypercharge field must 
be substituted and the resulting fluctuation projected onto the physical electromagnetic field. Similarly, 
this analysis will apply to any £7(1) fields in the low energy theory of inflation whatever their microscopic 
origin, although not to any non-Abelian fields. 
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is to study the leading departures from gaussianity, as measured by the bispectrum 
and trispectrum. Non-gaussian statistics have received considerable attention over the 
last several years. Indeed, recent experimental results hint that it may be possible to 
discriminate among different microphysical models using this technique [9j [10]. Such 
studies fall into the general framework of non-linear perturbation theory. However, 
departures of the observable adiabatic fluctuation from exact gaussianity are not the 
only effect we can hope to probe using this tool. In a quantum mechanical world which 
includes gravity we must expect the adiabatic fluctuation to couple to all other degrees 
of freedom which were light enough to be excited during the inflationary era. This point 
has been emphasized in recent work by Weinberg [TT]. If inflation does more for us 
than merely generate the observed density fluctuation — and if we are lucky — it may be 
possible to probe whatever other processes were in progress during inflation, provided 
they were sufficiently dramatic to leave traces in the inflationary observables. The 
analysis given in this paper can be thought of as an exploratory step in this direction. 

There are other motivations. If we are to have confidence in our predictions, it 
is necessary to maintain control over the theoretical tools which we use. This will 
become of increasing significance as the data improve and we pass from a qualitative to a 
quantitative description of the earliest times, excluding some models as possible theories 
of the early universe while promoting others as a better match for observation. The key 
criterion here is stability of the perturbative series which is used to extract observables 
from the Lagrangian, a question which has already attracted attention in the literature 
[T2| [13]. Several potential sources of instability exist. Increasing orders of perturbation 
theory are typically suppressed by powers of the ratio (H/Mp) 2 « 10~ 9 , where H is 
the energy scale of inflation and Mp is the Planck mass. However, perturbation theory 
can in principle contain instabilities which scale like a positive power of the scale factor 
a(t) fa exp(Ht) during inflation. 

If such "fast" divergences are present then they will rapidly overwhelm any powers 
of (if/Mp) 2 and render perturbation theory unstable after a short time — generally too 
short to be of any use in extracting predictions for relics of the early universe which are 
visible at the present day. Alternatively there may be large corrections which come from 
a sensitivity to physics in the ultra-violet, or from some other hierarchy which exists in 
the theory. The first possibility was studied by Armendariz-Picon et al. [15] . whereas 
an example of the latter, which was studied by Leblond & Shandera [T2] , is the relative 
hierarchy c<7 2 ^ 1 between the speed of sound and the speed of light. Whatever the 
source of large hierarchies which overwhelms the smallness of (H/M-p) 2 , it is important 
to emphasize that an instability in perturbation theory does not necessarily imply that 
anything untoward is taking place. It may simply mean that we need to find a better 
description of the process in question. 

The possibility of fast instabilities in perturbation theory was considered by 
Weinberg [HI [15] (see also Chaicherdsakul [16J), who gave a criterion according to 
which it is possible to decide whether such instabilities are forbidden. Even where this 
is the case, it does not necessarily follow that perturbation theory is convergent because 
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Weinberg's theorem does not exclude the possibility of much slower divergences: for 
any fluctuation which is outside the horizon, these divergences scale with the number 
of e-folds of expansion since the time of horizon exit. 

Whether divergences are fast or slow, however, the interpretation is the same. 
When we expand an expectation value of some operators as a series in a loop-counting 
parameter or the slow-roll parameters, we are developing a series expansion based on 
the cut-off associated with the theory. The role of the cut-off is played by the time at 
which we wish to evaluate the expectation value — which is usually chosen to be at the 
end of inflation, or some similar time where we wish to use the expectation value as 
an initial condition for classical cosmological perturbation theory in the later universe. 
The behaviour of any expectation value as a function of this cut-off is merely its time 
dependence. The problem arises because truncating the series at any finite order gives 
the appearance of divergences. If we compute an answer which is superficially divergent 
in this way, then we must find some other method to compute the time dependence of 
the expectation value before growing secular terms take perturbation theory out of our 
control. 

This point of view leads to a technique of computation in which we can separate 
calculations into a quantum initial condition [IT], for which we need all the machinery 
of the so-called in-in (or Schwinger-Keldysh) formalism, and a subsequent classical 
evolution for which to a good approximation we need only the classical, homogeneous 
evolution equations. The details of this approach have been developed by many authors 
[T8t [TFJ [19j Eni EED, EI]. There is a possible difficulty if we allow the inflating volume 
to become too large, because we may then encounter a source of quantum divergences 
which could invalidate the use of classical evolution equations even after horizon exit 
[22] 1231 EH E3 ESI E3 EH], but provided we work within some patch of spacetime not 
much larger than the size of the present horizon such effects are likely to be negligible. 

Most recent work on studying non-gaussianities from inflation has centred on the 
evolution subsequent to horizon exit [291 ESS EH [32], whereas the quantum initial 
condition has received comparatively less attention [331 EH EH EE]. There is a good 
reason for this imbalance: although there are known to be controlled examples where 
significant non-linearity can be generated outside the horizon [32l [36] EQ], it is very 
hard (with canonical kinetic terms) to construct a controlled calculation in which a 
significant effect arises from the initial condition. Indeed, the most useful tool for 
extracting predictions from the underlying quantum field theory — that is, the slow-roll 
expansion — generally has the effect of forcing correlations to be very small. A second 
interpretation of the analysis given in this paper is an example in which the initial 
condition is modified, by including a coupling to high-energy virtual quanta which 
belong to the electromagnetic field. As can be expected, it will not be possible to control 
the calculation in the regime where this modification dominates the initial condition. 
However, we will be able to obtain a bound on the characteristics of the interaction 
among scalar and electromagnetic quanta which guarantees that the calculation is not 
taken beyond our control. 
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Throughout this paper, we use units in which h = c = 1 and the reduced Planck 
mass is set equal to unity, giving M P = {&ixG)~ 1 / 2 = 1, where G is Newton's gravitational 
constant. The background space time is de Sitter space with flat spatial slices and metric 

ds 2 = -dt 2 + a(t) 2 dx • dx, (1.1) 

in the (—,+,+,+) sign convention. Some formulae are more conveniently written in 
terms of conformal time, defined locally by the rule dt = a(t) drj and given explicitly 
by the quadrature rj = f dt'/a(t'). Spacetime indices are labelled with Latin indices 
{a, b,c, • • •}; purely spatial indices are labelled with indices k, ■ ■ •}; and indices in 
the space of scalar fields are given Greek labels {a, (3, 7, • • ■}. 

Purely spatial vectors such as x or k are written in bold face and an infix dot denotes 
index contraction with the flat background spatial metric, so that x k = ^ Xiki = Xiki, 
where the summation symbol will usually be omitted. Note that both indices are 
lowered. This convention is used to interpret exponentiation of any square spatial 
matrix 7^, giving the rule exp(j)ij = Y^=o{l n )ij l n - Repeated spacetime indices 
in complementary raised and lowered positions are summed using the full spacetime 
metric g ab according to the Einstein convention, as usual; this convention also applies 
to raised and lowered spatial indices with the substitution of the full spatial metric in 
contractions. 

The model used as an example in this paper is Einstein gravity coupled to some 
collection of light scalar fields 4> a and a single U(l) gauge field A a . The U(l) gauge 
field is taken to have a kinetic term of the form \((j))F ab F ab , where F ab = d a A b — d b A a 
is the gauge-invariant field strength. The coupling, A(0), is determined by the vacuum 
expectation values of some or all of the light scalars. This model is studied in £j2j where 
the interactions among the scalar, tensor and gauge field perturbations are derived. In 
§3] the magnetogenesis mechanism [37 1 138 1 161 139] is briefly reviewed, following an analysis 
by Bamba & Sasaki. In §§1H5] I compute the leading loop correction — for the spectrum 
and bispectrum — which comes from scalar fluctuations mixing with virtual quanta of 
the gauge field. The paper concludes with a discussion in £j6j An auxiliary calculation 
of a simple pure scalar loop correction is given in Appendix A to aid comparison of the 
methods used in the present paper with those of other authors. 



2. Scalar— magnetic couplings in the inflationary Lagrangian 

Our starting point is Einstein gravity coupled to a collection of scalar fields a , some of 
which are responsible for inflation, with the addition of a U(l) gauge field A a . 

2.1. The gauge-fixed Einstein- scalar-vector action 

As discussed in the normalization of the gauge field is taken to be controlled by some 
non-canonical coupling A(0) which is determined by the vacuum expectation values of 
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some subset of the scalar fields. The action is 

S = i J d 4 x (r - V a Q V a Q - 2V - ^\((P)F ab F al }j , (2.1) 

where the £7(1) field strength is defined by F ab = d a A b — d b A a and R is the spacetime 
Ricci scalar. The potential V((j>) is any reasonably smooth function which will generically 
depend on all the <p a and is arbitrary except that in order for the analysis which follows 
to apply it must support an epoch of inflation, at least for some range of values of the 
scalar fields. 

The background spacetime is taken to be homogeneous and isotropic, so the gauge 
field has no expectation value, up to configurations which are pure gauge. One can 
therefore assume that A a is a perturbation, which will generically couple to the scalar 
and gravitational degrees of freedom in Eq. (12.11) . To study this system of coupled 
perturbations it is especially convenient to write the spacetime metric (including the 
effect of gravitational fluctuations) in the so-called Arnowitt-Deser-Misner (ADM) 
form, 

ds 2 = -N 2 dt 2 + h ij (dx i + N* dt)(dx j + N j dt), (2.2) 

where N (the "lapse function") and N l (the "shift vector") are not dynamical degrees 
of freedom, but instead are determined by constraint equations. The spatial metric hij 
contains propagating tensor modes, and depending on the gauge it may also contain 
propagating scalar modes. After gauge-fixing h^, and solving the constraints, N and 
iV* can be written in terms of the propagating degrees of freedom in and the <p a . 
This is a great simplification in concrete calculations. 

Consider first the gauge sector. The ADM decomposition is based on a so-called 
"3+1" split of spacetime into three spatial dimensions and one timelike dimension. Many 
formulae are simplified by making an analogous decomposition of the vector potential, 
writing A a = (p,uA where p is the timelike component, and uji is a spatial 3-vector. It 
will transpire that p is not a dynamical field, but is removed by a constraint associated 
with the gauge invariance of A a . Once p has been removed a further gauge fixing 
condition must be applied to u>i, which leaves the expected two physical polarizations 
of a massless gauge boson. 

When expressed in terms of p and u iy the part of the action involving the gauge 
field can be written 

SD^Jd 3 xdt NVh{-^h m U n f i3 f mn + j^hVdiudjU, | , (2.3) 

where we have defined a useful quantity BiU by the rule 

d t u = Ui- d lP + (2.4) 

(denoting a time derivative with respect to t by an overdot) and f\j is the spatially 
gauge-invariant 3-vector field strength, = diUj — djUi. The original gauge invariance 
associated with A a corresponding to U(l) gauge transformations was A a i— > A a + d a A 
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for some arbitrary spacetime-dependent function A. Under such a transformation, p 
and Mi undergo separate transformations determined by 



It follows that diU is gauge invariant, and therefore so is any action built out of and 



The ADM action for p and u>i is singular, because the lagrangian is degenerate 
along pure gauge directions. It therefore cannot be quantized as it stands, but must be 
put into a form suitable for quantum mechanical calculations by adding a term of the 
form s(c*f), where s is a so-called Slavnov operator, 



h is an auxiliary field, and 6, c* are a pair of anti-commuting ghost fields which 
nevertheless obey Bose-Einstein statistics. The function / is chosen to remove the 
degeneracy along gauge directions, but is otherwise essentially arbitrary provided that 
it is not invariant under gauge transformations. It is usually described as a gauge- 
fixing function. For the purposes of the present paper the most appropriate choice is 
an analogue of the Lorentz-Coulomb gauge, specified by / = cWjlj] It might have been 
thought that the covariant condition /' = d l u)i would be more appropriate in order 
to maintain manifest spatial covariance. Ultimately, however, we will be performing 
calculations in a version of the interaction picture in which it is simple to compute with 
/, but more complicated to compute with /'. It is immaterial whether we choose / or 
/', and since there are no serious consequences associated with losing manifest spatial 
covariance we will stick with /. (If desired, the spatially covariant action and Feynman 
rules can be obtained from those given here by the replacement <5 U h- > in the gauge- 
fixing terms written in Eqs. (12. 81) — ( 127251) below.) The total action is invariant under a 
quantized form of the original gauge symmetry, usually known as a BRST symmetry, 
which is generated by s. In the classical theory, the gauge condition is enforced by 
solving the constraint / = for one linear combination of the components of cjj. This 
removes one degree of freedom from the theory. 

After inserting the action into a path integral and appropriately integrating out the 
auxiliary field h, one finds that the ghosts b and c* decouple and contribute only to an 
irrelevant overall normalization. The result can be written as a path integral over the 
action 



p i-> p + A, 
tOi i— > oji + <9jA. 



(2.5) 
(2.6) 



diU alone. 







where the quantities B\ and _B_i are defined by 



Bl EE R - tfdi^drfa -2V- -tjVfijf, 



1 



SVS^diUjdmUn (2.9) 



B_i ee E ij E i:j - E 2 + 7i a 7i Q + Xh^diudjU. 



(2.10) 



f Recall that in the summation convention which is being used here, diU>i 
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In these formulae, R is the spatial Ricci curvature associated with h^; the quantities 
7r a = cf) a — NWj(f) a are a collection of momenta associated with the scalar fields; and 
Eij = hij/2— V(iiVj) is the gravitational momentum, where Vj is the covariant derivative 
compatible with h^. 

2.2. The constraint equations 

The physical degrees of freedom in the action (12.81) are a collection of propagating modes 
associated with the scalars, <f) a , together with modes arising from the components of hy. 
These are the fields whose time derivatives appear in the action. On the other hand, the 
fields N, N l and p do not appear in the action with any time derivatives, and therefore 
are associated with constraints. These constraints can be satisfied and the unwanted 
degrees of freedom eliminated by simply solving for N, N l and p in terms of the other 
fields in the system, and substituting the result back into the action. 

The constraints associated with iV and N l are well-known, and are modified here 
only by the presence of a component in the action associated with a gauge boson. The 
N constraint is 

R - h'^hcfOp,, -2V- ^//'"7/'"/, / ./;, i) . - ^6 mn d t u 3 d m u n 

- — {EijE i] — E 2 + n a n a + XWdiudju) = 0, (2.11) 
and the N l constraint is 

V* jl {E iS - M?) J = ^d 3 <p a - jjh ik diUf kj . (2.12) 
On the other hand, from integrating out p we obtain a very simple constraint 

V, (JjtfdiU^ = 0. (2.13) 

Consider first the N l constraint. We take the background field configuration to 
be spatially homogeneous and isotropic, so that the <ft a are functions of t only, with 
small perturbations 5<p a which satisfy the smallness condition \5<p a \ \(f> a \- If inflation 
has been ongoing for an exponentially large number of e-folds then this may be a poor 
approximation over the whole inflating volume, owing to back reaction effects which 
can cause large fluctuations in the scalar expectation values between widely separated 
regions. However, in any region of spacetime in the neighbourhood of exit from inflation 
this field theory is likely to be an accurate effective description. 

To parametrize the degrees of freedom associated with the spatial metric, we write 
hij = a 2 (t)(e' y )ij, where 7^ is a spatial 3x3 matrix whose components are taken to be of 
the same magnitude as b~(f) a . One then aims for a perturbative solution in 5<f) a and 7^, 
with the gauge field Ui also taken to be perturbative and of the same formal magnitude. 
This implies that the gauge field does not modify the background evolution of the scalar 
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fields. We define quantities a n , $ n and f3 n j by writing 

oo 

iV=l + ^« n (2.14) 

n=l 

oo 

N i = ^T l (d{& n + l3 nj ), (2.15) 

n=l 

where each of a n , d n and (3 n j is taken to contain exactly n powers of the perturbations, 
and djP n j = for all n. 

At 0(1), the N constraint (12.111) gives the Friedmann equation for the background, 

m 2 = l -j> a 4> a + v. (2.16) 

At 0(50), one obtains an equation for the scalar component of the shift vector, 
AH ■ 

—d 2 § 1 = -2VJ(j) a - 2<p a 5<p a + 2aA-6H 2 + 0°0A (2.17) 
a 2 

To obtain equations for a± and the vector component fiij one must return to the shift 
constraint, Eq. (12.121) . One finds /3y = and 

ai = (2.18) 

For the purpose of computing interactions between the S(f> a , •jij and u>i to third order, it 
turns out to be unnecessary to obtain any terms in the lapse or shift which are of higher 
order in the perturbations than Eqs. (I2.17l) -( [2.18l) . Although such terms are present in 
a naive expansion of the action, they are removed by the constraints (12. lip and (12.121) 
[35] . Eqs. (I2.17I) - (I2.18I) . together with the constraint /3%j = 0, are exactly the solutions 
for the lapse and shift which were found by Maldacena in a theory with no gauge boson. 
Indeed, at first order (although not above), neither the gauge boson u>i or any possible 
tensor modes 7^ contribute to iV or N l . 

To solve the p constraint one makes an analogous expansion of p in powers 
of the perturbations, writing p = Y^=i Pn with p n containing exactly n powers of 
perturbations. At 0(50) this gives 

d 2 Pl = diUi. (2.19) 
2.3. The Gaussian theory: Quadratic terms 

We are now in a position to study the evolution of the fluctuations 5<ft a , 7^ and 
cjj, together with their interactions. Any terms in the action which are linear in 
the perturbations vanish as a consequence of the background equations of motion. 
The leading non-trivial terms are therefore quadratic. Any theory defined by purely 
quadratic terms is free and gives rise to gaussian statistics. Therefore, the leading 
correction to gaussian statistics comes from interaction terms at cubic order or 
above. The details of these interactions can be calculated by treating them as small 
perturbations to the quadratic pieces, which we suppose still supply the dominant 
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evolution. This formulation is equivalent to the interaction picture in the canonical 
formalism. 

At quadratic order, the perturbations 5(p a , 7^ and Ui decouple and do not 
communicate with each other. The action therefore breaks into a sum of terms for 
each fluctuation which can be treated separately. 

Consider first the fluctuations in the light scalar fields S(f> a , which are described by 
the action 



5 2 3 - / d 6 x dt a' 



AH 



<9 2 tfi - 2V a 



20° 



aA-6H 2 + 



(2.20) 



In comparison the tensor modes 7^ obey a very simple action, containing only the 
minimal kinetic term 



So D 



d xdt a < jijjy ^dkjijdk'jij > . 



The action controlling the evolution of the gauge field satisfies 
S 2 ll \ I d^xdt a \\(uJi - dipi){ui - dip x ) - -^fijfij - -^diOJidjU 



(2.21) 
(2.22) 



2.4- Interactions: Cubic terms 



The leading interactions among the 7^ and a single scalar degree of freedom were 
obtained by Maldacena, and are not affected by the presence of a gauge field. However, in 
addition to the terms found by Maldacena there are now cubic interactions which involve 
gauge bosons together with the 5(p a or 7^. Because the spacetime gauge connexion 
A a appeared in the original action Eq. (12.11) quadratically, all these cubic interactions 
involve two gauge bosons and only a single other field. There is therefore a term which 
describes the interaction of two gauge fields with a scalar, 



S,dI [ d 3 zd*a 3 -f /V " n 



2A 

[uitbi - 2u i d i p 1 + dipidipi) + —ujif ij d j 'd 1 



\Jct) a 



2a 4 fok 



+ Cti 



* ifijfij ~ -rjdi^A^j - -4 (wA - 2u i d i p 1 + dipxdipx) 



>2.23) 



2a* JlJJlJ £a 4 

There is also an interaction between two gauge fields and a single 7^, 

5*3 5 ^ J d^xdt a 3 1 - (a;^ - 2u i d j p l + dip x djp x ) + ^"/ijfmifmj^ (2-24) 

In principle there is also a self-interaction between three gauge bosons, described by the 
term 



S 3 D i / iV.cdt <r 



- T UJ i diP2 

a z 



(2.25) 



We will see in £|3] below that it is possible to make a choice of gauge in which this 
interaction only involves the unphysical polarization of u>i. Therefore, in this gauge it 
decouples from all physical amplitudes, although it may remain present in more general 
gauges. 
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3. Magnetogenesis 

In this section the process of magnetogenesis is briefly reviewed, with the aim of 
establishing the relevant notation and formulae which will be required for a calculation 
of loop corrections in §§SH5j Our starting point is the assumption that the magnetic field 
arises from the coupling XF ab F a b between the gauge field and the scalar sector. This 
calculation has been given in some generality by Bamba & Sasaki [ID] , whose methods 
we follow. (See also Refs. pEH M, S3].) 

When calculating processes involving exchange of virtual gauge bosons in 
Minkowski space it is usually a useful calculational check to leave the constant £ 
arbitrary. Indeed, since physical quantities do not depend on £ it must cancel out in 
any correct computation. Unfortunately, when attempting to compute the gauge field 
propagator in a time- dependent background, the presence of the £ term is an obstruction 
to solving the propagator equation. To simplify the process, it is helpful to take the 
limit £ — > 0. For a gauge-fixing functional / this directly enforces the constraint / = 
in order that the action remain non-singular. Hence, instead of integrating over the 
three components of u)i, we must make the decomposition 

"(CM) =J2f ^<(kK(M)e ikoc , (3.1) 

where the e a are so-called polarization vectors, labelled by a two-valued discrete index 
a. The path integral should include only the two physical polarizations u±. The 
polarization vectors are chosen in such a way that k • e CT (k) = for each a, and are 
normalized so that 

e CT (k) .e a '(k)* = 6 mT> (3.2) 
£ e f(k) e j(k)* = ^-^. (3.3) 

<re± 

After carrying out this reduction, note that the apparent three-boson interaction (12.251) 
couples only to the unphysical degree of freedom in It can formally be removed 
after integrating by parts. It follows that this interaction is not physical. Note also that 
Pi = in this gauge, so all p terms drop out of the action for uji up to cubic order. 

With these choices, the propagator for the gauge field is obtained by inverting the 
differential operator which appears under the integral in Eq. f)2.22p . Suppose we write 

(w (T (k 1 ,t)uv(k 2 ,f')) = (27r)5(ki +k 2 )<WG fe (M')- (3.4) 
The equation which determines G k is 

G k + (h + ^ j G k + ^G k = -±S(t - 0, (3-5) 

together with the boundary condition that in the limit k/aH — > oo, where the 
fluctuation corresponding to this wavenumber cannot feel the curvature of spacetime, 
G should approach the corresponding mode function from Minkowski space. This 
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boundary condition plays an important role in the calculation. It corresponds to the 
stipulation that we begin with conventional Minkowski space quantum field theory in 
the ultra-violet, and then attempt to determine what this implies for fluctuations deep 
in the infra-red. Although it is possible to make more general choices of field theory in 
the ultra-violet, the assumption of flat space field theory is minimal and any admixture 
of different ultra-violet physics is usually subdominant to the Minkowski result. When 
we come to define what we mean by loop integrals in de Sitter space, which are also 
part of the specification of the ultra-violet behaviour of the theory, it will be necessary 
to take care that our definition does not destroy the property that we begin with flat 
space field theory at very high energies. 

Eq. (I3.5P cannot be solved exactly for a general choice of A (</>). Instead, one can 
obtain a solution of Wentzel-Kramers-Brillouin (WKB) type which is valid inside the 
horizon, and can be matched onto a long wavelength solution which is valid outside the 
horizon. This is equivalent to using the flat space boundary condition deep inside the 
horizon to determine the size of fluctuations at horizon exit, and then using this as an 
initial condition for a classical superhorizon calculation. 

The solutions can be written most simply in terms of the conformal time coordinate 
f]. In this variable, the relevant WKB solution is 



where V denotes evaluation at the time the mode with wavenumber k exited the horizon. 
The k dependence gives the spectrum a steep blue tilt, which means that if it remains 
unprocessed by new physics in the superhorizon regime it must have an essentially 
negligible magnitude on observable scales. 

It is sometimes said that a canonically normalized vector field does not receive 
a perturbation from inflation, based on the observation that positive and negative 
frequencies of the gauge field are not mixed as the universe expands [HJ 05] . It follows 
that if asymptotic in and out vacua can be defined there is no particle creation in 
the transition between the in- and crai-vacuum. In the computation of inflationary 
observables, however, there is usually no natural out region where a notion of particles 
make sense, nor any need to invoke such a region. Instead, we are interested primarily 
in whether expectation values of operators behave coherently over many e-folds of 
expansion, and in this sense a canonically normalized gauge field receives a fluctuation in 
exactly the same way as any light bosonic field. As we have already observed, however, 
the fluctuation which is imprinted in the precisely canonical case is very blue and entirely 
negligible on cosmologically relevant scales. One can think of this as a consequence of 
the fact that canonically coupled vector fluctuations must redshift like radiation, giving 
an extremely strong suppression for scales which exited the horizon early and have been 




7] < 7] 
7]' < 7] 



(3.6) 



It follows that the power spectrum of each polarization at horizon exit satisfies 




(3.7) 
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redshifting for longer. Once modes re-enter the horizon they begin to oscillate and 
expectation values of their associated operators lose their coherence, which is consistent 
with Parker's observation that there is no asymptotic particle creation in this model 

The situation changes in the presence of the non-canonical coupling A(0). In this 
case Eq. (13. 7p suggests that if A* < 1 the power in fluctuations of the gauge field 
has been amplified in comparison with a vanilla model where A = 1 for all time. 
However, this is misleading. Since the energy-momentum tensor associated with the 
gauge field is proportional to A, when we compute the energy density stored in Ui at 
horizon exit we obtain an answer which is independent of A*. It follows that there is 
a subtle distinction between this method of generating perturbations and the familiar 
case of scalar perturbations. In the simplest model of scalar perturbations, we generate 
fluctuations of the correct magnitude at horizon exit which are then conserved until 
horizon re-entry. In the case of Ui we do not make the physical fluctuations at horizon 
exit any larger whether or not we include a coupling A. It is evolution of A after horizon 
crossing which amplifies fluctuations in the gauge field and prevents them redshifting to 
zero; in this respect, the mechanism is similar to the curvaton example for purely scalar 
perturbations. We should therefore expect physical quantities to involve the ratio, A 2 / Ai, 
of A at different times 771 and 772, which can be made large only if A evolves sufficiently 
strongly that there is a large hierarchy between its values at these times. In the limit 
of an observation made instantaneously at a moment in time, we can conclude that the 
physical effect must be proportional to A'/A, or a higher derivative, where ' denotes 
differentiation with respect to the conformal time. The power spectrum and bispectrum 
measured at horizon exit are examples of observations made instantaneously, and we 
shall see in §§1H5] below that they come proportional to powers of A'/A. An alternative 
means of breaking conformal invariance was considered in Ref . |16] . 

Once a mode has left the horizon, its evolution is governed by Eq. (13.51) in the 
limit k/aH — > 0. Discarding powers of gradients, the homogeneous solution for each 
polarization mode takes the form 



which depends on the value of u and its derivative u' (where a prime ' denotes 
differentiation with respect to conformal time) on any initial hypersurface, 77 = 77*, 
provided that all relevant modes have left the horizon at that time. For a mode 
corresponding to a single wavenumber k we can take this hypersurface to be the time of 
horizon exit, and the power spectra of 10 and u' can be extracted from Eq. (13.61) . If A is 
increasing or decreasing less fast than 77 then the integral in Eq. ( 13. 81) converges and each 
polarization is constant outside the horizon up to terms which decay exponentially fast 
in cosmic time. On the other hand, if A decreases faster than 77 then the integral diverges 
and each polarization grows rapidly. These possibilities correspond to a decreasing gauge 




(3.8) 



f I would like to thank David Lyth for helpful correspondence on this question. 
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coupling (or one which increases less fast than a -1 ), or an increasing gauge coupling, 
respectively. 

At any time r\ at which Eq. ( 13. 8ft applies, the proper electric and magnetic energy 
densities were computed by Bamba & Sasaki and are given by [HJJ 



B^^e M v) = ^(yx^, (3.9) 

where is the alternating tensor in three dimensions and the normalization has been 
chosen so that the electromagnetic energy density is given by pem = (B 2 + E 2 )/2, as 
usual. It follows from Eqs. (I3.9I) - (I3.10I) . Eq. (13.81) and Eq. (13.61) that the power spectrum 
of the proper magnetic energy density is 

1 A / k" 



An 2 A* V a 



2 XJ J, AO 



and the power spectrum of the proper electric energy density is 



„ 1 K k 2 
ri 



1 A' 2 



iib- 

2 A, 



(3.11) 



(3.12) 



4vr 2 A a 4 

Bamba |47J (see also Giovannini [12] and Martin & Yokoyama [IS]) observed that if A 
is increasing (or decreasing less fast than a -1 ), then the integral is dominated by its 
lower limit and one finds that Ve/Vb ~ (A*/A) — > as rj — > 0. On the other hand, if 
A is decreasing sufficiently fast to cause the integral to diverge then it is dominated by 
its upper limit, giving instead Ve/Vb ~ (aH/k), which grows exponentially with the 
number of e-folds since horizon exit. One can conclude that in the first case one has a 
predominantly magnetic field at late times, with only an exponentially small admixture 
of electric field, whereas in the second case the situation is reversed. 

Let us focus on the case where a magnetic field is synthesized at late times, 
reserving the opposite case for the discussion in £j6l During radiation domination the 
electromagnetic field redshifts like the dominant constituent of the universe and therefore 
its relative density is conserved. If we assume prompt reheating after inflation, then the 
root mean square fluctuation in on a scale corresponding to comoving wavenumber 
k has magnitude 



IBU-^-J (3.13) 

in Planck units. The proper wavenumber today on cluster scales is of order k/a ~ e -140 , 
and to seed a galactic dynamo it may be sufficient to produce fluctuations with 
magnitude [49J |B| 

rms ~ 6 T|ij Therefore to obtain a cosmologically interesting 
magnetic field, we require roughly A ~ e 200 A* [IE]. Although the mechanism of 
magnetogenesis is quite insensitive to the dependence of A on the scalar fields which 

| Our conventions for magnetic field strengths are measured in tesla, where IT ~ e -120 in Planck 
units. 
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determine its value, a large class of models which invoke couplings of this form can 
be written in the "dilaton-like" form \(<p) = exp((f> / M^) [50], where is some 
characteristic mass scale. Such couplings have also been invoked in the context of models 
of dark energy [51], where they may be subject to additional constraints [52, 53] . In 
order for a coupling of this form to yield the requisite hierarchy, must be chosen to 
satisfy « A0/2OO where A</> is the excursion of the field <fi between horizon crossing 
and the end of inflation. If is a field driving a stage of chaotic inflation it can be 
supposed to move a distance perhaps of order 10 in fundamental units, in which case 

~ i (3.14) 

We will adopt this value of as the canonical one for magnetogenesis, although in 
practice M$ will vary from model to model, and may be closer to the Planck scale. 

A variety of bounds are known on the energy density which can be present in 
magnetic fields at early times, either from direct detection [54] or indirect effects [55] . 
These limits typically arise from constraints at the short wavelength end of the spectrum 
and imply that the extremely blue raw spectrum, Eq. (13.71) . must be processed into an 
approximately scale-invariant form. Exact scale invariance occurs for = v2e/4, 
which gives M^ 1 ~ lojjj and it follows that for close to Eq. (13.141) approximate scale 
invariance will apply. 



4. Electric loop corrections to the scalar spectrum 

In this section we return to the interaction between gauge bosons and the other 
fluctuation modes which are relevant during inflation. If we wish to make predictions 
for the anisotropy seen in the cosmic microwave background (CMB) then we are 
principally concerned with the power in scalar perturbations around the time of last 
scattering, because these dominate the density fluctuation in the primordial plasma. 
Such fluctuations are connected to observation by making predictions for the properties 
of the comoving curvature perturbation, TZ, which on superhorizon scales is equivalent to 
the curvature perturbation on uniform density hyper surf aces, (. In a model with many 
degrees of freedom there is a considerable simplification afforded by using (, which can 
be computed using the so-called non-linear 5N formalism. 

As an initial condition, calculations using the 5N formalism require predictions for 
the correlations among the 5<f) a around the time of horizon crossing. In this section 
we will compute a prediction for the two-point correlation of the 50°, taking into 
account the leading loop correction from exchange of virtual gauge bosons. This is 
presumably sufficient to make predictions for the power spectrum of the CMB, although 
it will transpire that there may be UV-divergent terms associated with the gauge 
transformation between ( and the 5<p a which are not captured by the 5N formula. 
Therefore the final answer can be treated as an order of magnitude estimate only. 

§ Compare, for example, with Eq. (42) of Ref. [38] . 
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In the following section we will compute the analogous correction for the three-point 
correlation function, which is necessary if we wish to study higher-order statistics. 

4-1- Dominant contributions to the interaction vertex 

Consider the vertex for interaction of two gauge bosons with a scalar particle, given by 
Eq. (12.231) in the limit £ — > with Ui replaced by Eq. (13.11) . If we are only computing 
around the time of horizon crossing, then we can obtain a good approximation by 
truncating all quantities to leading order in the slow-roll expansion. In the interaction 
vertex, this leaves us with terms of the form 

S3 12 2 J d 3 xdt a X a 5(j) a yliiUi -{diUJjdiUj — diUJjdjUi) j . (4.1) 

The first term in Eq. (14.11) involves an interaction with uj 2 , which according to 
Eq. (I3.10p can be thought of as a measure of the electric field intensity. The second 
term involves interactions with terms of the form (du>) 2 , which according to Eq. (13.91 ) 
are a measure of the magnetic field intensity. We can likewise imagine interactions which 
are dominated by the first or second term to represent interaction with the electric or 
magnetic field, respectively. This distinction is useful because if A is carrying a strong 
time dependence at horizon exit we expect the ambient electric field to be enhanced in 
comparison with the magnetic field, which depends only on spatial gradients. Indeed, 
this is true irrespective of whether A is increasing or decreasing, provided that it has a 
strong time dependence in either direction, because the properties of the fluctuations at 
horizon crossing do not determine whether the final configuration will be an electric or 
magnetic field. The initial condition involves a strong electric field in either case, but 
the final character of the field is only determined by the long-term evolution of A after 
horizon crossing. 

To see this in detail, it is simplest to use Eqs. (13.11) and (13.61) to rewrite Eq. (14.11) 
as an effective vertex which takes the form 
f , f d 3 A;i d 3 A; 2 d 3 fcs . 

53 " J ( 27r )9 ( 27 <K k l + k 2 + k 3) 



A 

r 

where the vertex function ct % L is defined by 



x '-jSrO^i, /7)4 ± e?(k 2 )e;(k 3 H(k 2 , V )u b (k 3 , r,), (4.2) 



a% l3 = 5 %3 jk 2 • k 3 + (u 2 k 2 + ^Mj (iL 3 k 3 + J - k 2j k 3l . (4.3) 

In this equation, Q is an abbreviation for the dimensionless combination 

a = (4.4) 

2 A H V ; 

and the ± symbols are chosen according the details of time ordering and the assignment 
of '+' or '— ' vertices (to be described in §4.21 below) among the gauge fields which 
participate in the vertex, and are fixed by the structure of the diagram in which Eq. (14.21) 
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is inserted. (In Eq. ( 14.31) . we have temporarily abandoned our summation convention — 
for this equation only — in the interests of notational clarity: this choice of a*^' should 
be inserted directly in Eq. (14.21) without concern for the position of the spatial indices i 
and j.) The parameter Q measures the hierarchy between the wavefunction of the gauge 
field and its time derivative, and arises from the interaction with the electric field. In 
the limit ^> \kr)\, the total interaction is dominated by this electric piece and the 
details of the time ordering become irrelevant. This is true whenever the rate of change 
of the coupling A with the fields is tuned to satisfy 

>e" JV *. (4.5) 

In Eq. (14.51) . '*' denotes evaluation at the time when the mode with wavenumber 
k exited the horizon, or more precisely a small but non-zero number of e-folds iV* ~ 1 
afterwards. This time should be chosen so that the fluctuations in scalar modes are 
close to their asymptotic superhorizon values, and the canonical commutation relation 
[5(f), S(j)} ~ e~ N * allows the 5<ft a to be treated as approximately classical quantities, but it 
should not be so late after horizon exit that appreciable evolution might have occurred, 
which would spoil the accuracy of the slow-roll approximation. Whenever Eq. (14.51) 
applies we can ignore the purely magnetic part of the interaction. It is important to 
note, however, that the electric interaction can only become dominant if at least one 
of the scalars which couple to u>i is rolling on cosmological timescales, so that <fi a /H is 
not totally negligible for some <p a . This is not really a restriction if one wishes to use 
this interaction for the purposes of magnetogenesis, because one is then dependent on 
A developing a large hierarchy between its value at horizon exit and its value at some 
much later time, such as the end of inflation, and this can occur only if some of the 
scalar fields are in motion. 

Alternatively, the scalar fields whose vacuum expectation values determine the 
magnitude of A may not be rolling during inflation, or the dependence of A on these fields 
may be so weak that (lnA) a is never large enough for Eq. (14. 5p to apply. In this limit 
the interaction is still interesting, but it is magnetically dominated and the momentum 
integral which describes the loop is somewhat more complicated to compute. For the 
remainder of this paper, we focus on the electric part of the interaction only. 



4-2. The in-in formalism 

The appropriate formalism in which to compute expectation values of cosmological 
fluctuations is the so-called in-in formalism introduced by Schwinger. In the 
cosmological case we wish to compute expectation values of the form (0|O|0) for some 
collection of operators O, in the state |0) which following the discussion of ultra- 
violet physics given in §|3] should be taken to be the Minkowski vacuum deep inside 
the horizon. We know from scattering calculations in Minkowski space that one can 
compute (out|0|in) using a conventional Feynman path integral for any 'in-state' |in) 
and 'out-state' |out). It follows that after integrating over all possible out-states we can 
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compute (in | O | in) using two path integrals, which gives the so-called Schwinger-Keldysh 
path integral formula 



where the label the elementary fields in the theory and O is taken to be built out of 
either '+' or ' — ' fields but not both. The integral is defined by prescribing that only 
those + and — field configurations which begin in the appropriate vacuum |0) and end 
with coincident values at some late time are to be included in the integration. The 
precise choice of this late time is immaterial, provided it is chosen to be later than the 
time of evaluation of any field which appears in O. We will conventionally choose O to 
be composed only from + fields and carry the integral from past infinity to the time of 
observation, 77,. 

When coupling to gravity is taken into account, there can be a subtlety in the 
construction of the path integrals which appear in Eq. (14.61) . From a given formula 
for the lagrangian, one would ordinarily obtain the canonical momenta and construct 
the hamiltonian. Integrating over the coordinates and canonical momenta, with time 
evolution supplied by the hamiltonian, one arrives at the standard path integral. If 
the lagrangian is not quadratic in the momenta, however, then one cannot explicitly 
integrate them out |56j. Instead, one must include their degrees of freedom in the path 
integral as "derivative ghosts" which compensate for the fact that the fields are not 
canonically normalized. They are associated with interactions which are cubic or higher 
in derivatives such as 8(f). In our example, the only such interactions are associated with 
scalar fluctuations and for this reason we ignore derivative ghosts. Instead, they should 
instead included with scalar loop corrections, which presumably do not lead to large 



The doubling of degrees of freedom in Eq. (14.61) implies that when we expand the 
path integral into diagrams we encounter extra graphs, which mix vertices constructed 
from + and — fields. These vertices collectively ensure that all expectation values 
are real. If we apply Eq. (14. 61) to the calculation of the scalar two-point function 
(50 a (ki)50 /3 (k2))*, one obtains the diagrams shown in Fig. [U 

The diagrams in Fig. [1] divide into two pairs of complex conjugates, corresponding 
to the (+, +) and (— , — ) diagrams in one pair and the (+, — ) and (— , +) diagrams in 
another. Consider first the (+, +) diagram. This makes a contribution to the scalar 
two-point function of the formpj 



f In writing these and all subsequent expressions, I have chosen to nest the integrals which result from 
mixed (+, — ) and (— ,+)-type diagrams. Alternatively, it is possible to factorize these contributions, 
obtaining expressions which are manifestly non-singular for all momenta ki [591 1601 161j . Whichever 
method is chosen, the answer is always the same. I would like to thank Peter Adshead, Richard 
Easther, Eugene Lim, Martin Sloth and Filippo Vernizzi for correspondence on this issue. 




(4.6) 



effects [561 EH EE]. 
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Figure 1. Schwinger-formalism diagrams for the loop correction to the scalar two- 
point function which arise from mixing with gauge bosons. Straight lines indicate 
scalar quanta, which only appear on the external legs. The interior loop, composed of 
wavy lines, indicates mixing with virtual quanta borrowed from the ambient electric 
and magnetic fields. The fields associated with external legs are always of + type, 
whereas the vertices can be of + or — type. The (+, +) and (— , — ) diagrams form one 
complex conjugate pair, and the (+, — ) and (— , +) diagrams form another. 
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(4- 



To make sense of this for real (3\ and /3 2 , the contours of integration for rj and r must 
be deformed so that e lt3lTJ and e 1/32T are decaying for large \rj\ and |r|, respectively. 
Indeed, this contour prescription follows automatically from the choice of vacuum |0) in 
Eq. (14.61) . In addition, P2(q, r ) is a polynomial in the vector momenta which depends 
on the transverse structure of the gauge field propagator and the detailed momentum 
dependence at the vertex. It satisfies 

P 2 (q,r) =gV + (q-r) 2 . (4.9) 

Eq. (14. 7p has been written in a form which emphasizes the symmetry between k x and k 2 . 
In order to simplify it further, one can integrate out either q or r, which leaves behind 
a single momentum-conservation delta-function, S(ki + k 2 ), together with an integral 
over the remaining momentum. Since local field theories usually exhibit bad behaviour 
at high energies, we may expect this integral to receive a significant contribution from 
the region where the 3-momentum which circulates in the loop becomes large. Indeed, 
neglecting any powers of momentum which are introduced by the time integrals in 
Eq. (14. 8p . it easy to see that the loop diverges at least as fast as J dq in the ultra-violet, 
whereas it converges at least as fast as J qdq in the infra-red. The time integral has 
dimension q 2 , so it can only make this divergence worse at high energy, and improve 
convergence at low energy. To make sense of such an integral it must first be regularized, 
removing the contribution of arbitrarily high-energy modes, and made finite by applying 
a renormalization prescription. Once this is done, the ultra-violet sensitivity of the 
original integral implies that we can expect the loop to be dominated by contributions 
near some UV scale. In what follows, we use this to ignore effects associated with the 
infra-red cutoff. 
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How are we to choose a cut-off in an expanding, inffationary spacetime? In flat space 
quantum field theory we are familiar with the use of a variety of regulators. The method 
of the renormalization group tells us that these are all equivalent in the continuum limit, 
after the introduction of counter-terms and a renormalization prescription. The simplest 
choice is a hard cutoff, M, on the invariant momenta associated with particles which 
circulate within the loops. To apply this cutoff in practice one must Wick rotate loop 
integrals to Euclidean spacetime where M becomes a SO (4)-invariant cut-off on the 
Euclidean momentum. When we revert to Lorentzian signature, this procedure gives 
a 5*0(3, l)-invariant result and therefore preserves Lorentz invariance. As a result, the 
hard momentum cutoff in Minkowski space provides a popular means of estimating the 
magnitude of loops in an effective theory. The problem at hand is to find a way to 
perform similar estimates in de Sitter space. 

In any curved spacetime we do not have global Lorentz invariance, although 
according to the equivalence principle we must recover approximate local SO (3,1) 
invariance in a small neighbourhood of any point. On the basis of the discussion in §3] 
it follows that this restoration of Lorentz invariance in the ultra-violet corresponds to 
choosing conventional flat space quantum field theory deep inside the horizon. Whatever 
regulator we pick, it is necessary to be sure that it does not conflict with local Lorentz 
invariance in the ultra-violet, or we shall obtain nonsensical results. Recently, van 
der Meulen and Smit [62J have pointed out that that a momentum cut-off at some 
characteristic scale M should be taken to apply to proper momenta, rather than 
comoving momenta, since only the former have physical significance. However, such 
a cut-off violates local Lorentz invariance for any finite value of the cut-offl]] This is 
because the cut-off grows between adjacent spatial slices which are taken at later and 
later times. If we pick some particular point in spacetime and ask whether Lorentz 
invariance is restored in a local neighbourhood of that point we find that no matter how 
small a patch of spacetime we choose, it necessarily intersects a sheaf of spatial slices. 
The growth in the cut-off transverse to the sheaf picks out a preferred direction, and 
the properties of the theory do not become invariant under local spacetime rotations 
within the patch. It follows that Lorentz invariance is broken. 

Breakdown of Lorentz invariance in the ultra-violet would have dramatic 
consequences for the stability of the theory. As time increases, the cut-off grows. This 
means that between two adjacent spatial slices, one taken at a slightly later time than the 
other, new quanta are introduced to our universe [66]. These new quanta destabilize 
the theory, because although they are at very high energies when they are added to 
the description, they can scatter with each other to produce other quanta which are 
much softer. These scattered soft quanta, whose correlations are essentially random, 
impinge on superhorizon correlation functions and scramble their values. Thus, in the 
presence of the Lorentz-violating cut-off, correlation functions do not retain coherent 
values outside the horizon but instead behave as if they are coupled to a quantum noise 



| I would like to thank Andrew Tolley for emphasizing this property of the proper momentum cut-off. 
A similar observation has been made in several places in the literature; see especially Refs. [63 | 164 } [65 ] . 
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bath and are swamped exponentially quickly by quantum noise. In this cut-off theory 
there is no quantum-to-classical transition, and presumably there are no conserved 
quantities outside the horizon. On the other hand, if we find a way to preserve Lorentz 
invariance then none of these undesirable features become manifest. One can think of 
this stabilizing property of Lorentz invariance as analogous to the conservation of the 
velocity of the centre of energy of any isolated system. This is a consequence of invariance 
under Lorentz boosts and implies, for example, that quantum effects associated with 
virtual ultra-violet quanta do not disturb the motion of isolated particles. A form of 
this analogy was earlier used by Lyth [67J. 

Let us return to Eqs. (14. 7B - (14. 8 p and study the time and loop integrals in more 
detail. Along the way, we will encounter the features described in the previous 



paragraph. In Appendix A an analogous analysis is given for the scalar loop which 
arises from a V" interaction, which has previously appeared in the literature [621 125] . 
The (+, +) diagram described by Eq. (14.71) must be supplemented by a (+, — ) diagram, 
which gives a contribution corresponding to 
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(4.10) 



Once we have aggregated the contribution of Eq. (14.71) and Eq. (I4.10p we must remember 
to add in their complex conjugates in order to account for the (— , — ) and (— , +) 
diagrams which have not been written explicitly. 

After integration by parts in rj and r, the fundamental time integral given in 
Eq. (14.81) can be re-expressed in the form 
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(4.11) 



This equation shows primitive fast divergences but some of these will cancel and others 
are purely imaginary, which implies that they disappear when the contribution of all four 
diagrams is accounted for. Write k\ = k 2 = k and consider the deep ultra-violet region, 
where q ~ r ^> k. In this region, the fa terms are individually very large, approximately 
satisfying $ ~ ±2g ^> k, so that the integrals are almost all very small — except for 
the integral with (3\ + (3 2 in the exponent, since this combination is independent of gj§| 

§ Conversely, these integrals diverge like powers of In \qr)* \ near the infra-red cutoff, where q is a measure 
of the 3-momentum circulating in the loop. However, in this region such logarithms are suppressed 
by positive powers of q and therefore the loop integral will be well-behaved. Indeed, for the purposes 
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Collecting all necessary terms, expanding for small \kr]*\ and keeping only contributions 
which are relevant in this limit, the expectation value can be written in the form 

(5r(ki)5/(k 2 )), = (27r) 3 5(^k^^^-(^) Q%, (4.12) 

with the left-over time and momentum dependence consolidated into an integral of the 
form 

h = 4 J g 3 |k- q|3 P2(q ' k ~ q) ( fc2 ~ Hq + |k ~ q|)(iV * + ln2 - 7e - 1)}, (4.13) 

where 7e ~ 0.57722 is the Euler-Mascheroni constant and iV* measures by how many 
e-folds the mode with wavenumber k is outside the horizon at time 77*. 

Eq. (14.131) is quadratically divergent. Introducing a cut-off A (as yet unspecified) 
on the momentum circulating in the loop, this integral takes the form 

Rk Ah 2 Ifit 3 A 

h « -^A 2 (iV* + In 2 + 7E - 1) + —A + ^-ln -(N m + ln2 + 7e - 1) + • • • , (4.14) 

where '• • •' denotes terms which are subdominant in the limit A — > 00, and \x is an 
arbitrary scale, of the same dimensions as A, which has been introduced to make sense 
of the logarithm. It is clear that if we take A to correspond to a proper momentum 
cut-off, which would take the form A = a*A for some constant A, then we introduce 
fast divergences as a t -> 00. These divergences correspond to hard quanta which are 
redshifted into the effective theory as the universe expands but subsequently scatter 
to produce soft quanta and contaminate the spectrum, as described above Eq. (I4.10p . 
Indeed, since the accumulation of such quanta is presumably highly incoherent, it seems 
unlikely that one can ascribe any definite value to the spectrum after horizon crossing. 

The key lesson I wish to draw from this example is that such divergent effects 
are fictional. If we begin with a Lorentz invariant theory valid at high scales and 
integrate out modes above a proper momentum cut-off, then cancelling divergences 
would automatically appear in the coefficients of the resulting effective Lagrangian 
[68, 69J. For this reason, power law divergences such as those appearing in Eq. (14.141) 
are devoid of physical significance; only the coefficient of the logarithmic divergence can 
have meaning. If we begin with a Lorentz invariant effective low energy theory such as 
Eq. (12. ip it is not possible to see these cancellations taking place. Thus, taken literally, 
our analysis up to this point is not compatible with a Lorentz violating cut-off; instead, 
sensible answers can be obtained only by using a Lorentz-invariant regulator such as 
dimensional reduction rTJJ [15], [16], [58], [60] . However, we can equally well make use of 
our knowledge that the power law terms in Eq. (14.141) must ultimately cancel, leaving 
behind an unfixed finite term or threshold correction [70] [71]. 

The threshold correction can be determined by matching to a more complete theory 
which resolves the details of ultraviolet physics, or by specifying a renormalization 

of the present calculation, we are assuming that the loop integral is dominated by exchange of virtual 
quanta near the ultra-violet cutoff, so these integrals (and other similar integrals to be encountered in 
S}5] while studying the loop-corrected three-point function) actually play no role in the analysis. 
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prescription which allows us to make contact with measurement. However, unlike 
simple theories such as quantum electrodynamics it is difficult to find an appropriate 
renormalization prescription in cosmology. This is because it is not possible to directly 
measure the expectation values we have computed: they are only important as an 
initial condition for the purpose of computing the structure in the late universe which 
is visible to us. A similar problem afflicts calculations in quantum chromodynamics, 
where interactions among hadrons are handled by first studying the predominantly 
electromagnetic interactions among their constituent partonsjjj] In the case of QCD it is 
also impossible to observe the parton correlation functions directly. 

Without extra information, little can be concluded about the magnitude of any 
possible threshold correction. If it is large it would correspond to a non-negligible 
mass for the inflaton, and therefore reproduces the well-known ^-problem. Since it 
is necessary in any case to assume that some conspiracy of ultra-violet physics allows 
inflation to occur at all it is presumably a reasonable approximation to set the threshold 
correction to zero. We can then estimate a lower limit for the magnitude of the 
loop correction by supposing that the logarithm receives contributions up to the scale 
where renormalization begins to remove ultra-violet modes from the theory, leaving a 
contribution of the form ln(A'//x), where A' is the scale of new physics. The nature of 
the cutoff is now immaterial. Taking /i « k « a*H* and A' to be the proper scale of 
new physics, it follows that the loop corrected spectrum can be written 



where Po(k) = H 2 /2k 2 is the tree-level power spectrum, and Vo(k) = H 2 /4n 2 is its 
dimensionless equivalent. Note that since a flat metric is being assumed on field space, 
the placement of the a and f3 field indices is immaterial. Eq. (14.151) is the first principal 
result of this paper. It is interesting to observe that the logarithm hx(A'/H^) cannot 
be too large; although the precise value of during inflation is model dependent, in 
chaotic models it is reasonable to assume that H* ~ 10~ 5 in fundamental units. It 
follows that ln(A'/f/*) cannot be more than of order 1 - 10. 

The two-point function (<50 Q (ki)5<^(k2)) is one contribution to the power spectrum 
C, but in practice it would be accompanied by other contributions arising from non- 
11 Indeed, there is an interesting analogy between quantum chromodynamics and the calculation of 
cosmological expectation values using the SN formalism. In QCD, one calculates correlation functions 
among partons at high energies, where the QCD coupling is small and perturbation theory applies. One 
then chooses a "factorization scale" which determines the energy below which partons are summed 
up into hadrons according to certain "parton distribution functions." In doing so, it is possible to 
encounter large logarithmic singularities at soft or collinear momenta which enhance the phase space 
for interacting partons to dress themselves into jets, or which manifest as initial state radiation. One 
can find analogues for all these effects in cosmology: high energy corresponds to early times during 
inflation, where the slow-roll approximation applies and perturbation theory in slow-roll quantities 
makes sense; the time of horizon crossing plays the role of the factorization scale; the coefficients of the 
SN expansion correspond to the parton distribution functions; and logarithms such as iV* = ln|fe?7*| 
play the role of the singularities at soft or collinear momenta. The details of this analogy have been 
explored in more detail in Ref. [72] • 



P af, (k)=Po*(k){s a < } + ~V *tt 




(4.15) 
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linear terms in the gauge transformation between £ and the 5(j) a . These may themselves 
carry ultra-violet divergences which should be accounted for in an accurate calculation. 
However, there is no reason to believe that these contributions would be any larger than 
Eq. (14.151) . and we can therefore suppose that this expression suffices for the purpose of 
obtaining order-of-magnitude estimates. 

How large is the loop correction? Specializing for simplicity to the case of a single 
field, and adopting the parametrization \(<p) = exp(4> / M^) , it follows that this loop 
correction does not overwhelm the tree-level provided M~ l < 120. This compares 
favourably with the value M^ 1 < 20, given in Eq. (I3.14p . which was suggested as 
appropriate for magnetogenesis (see ^3} although M^ 1 varies from model to model, it is 
unlikely to be as large as 120). Note that although these values may seem closer than 
desirable, the tuning here is in an exponential. A mass scale M^ 1 = 120 corresponds 
to a hierarchy A/ A* = e 1200 which is enormously larger than is required or desirable 
to produce a primordial magnetic field: the large energy density of an electromagnetic 
field amplified by such an enormous factor would swamp any other constituents of the 
universe and lead to an entirely unacceptable late time phenomenology. 



5. Electric loop corrections to the scalar bispectrum 



Eq. (I4.15P and the associated bound on A or which guarantees that the loop 
correction does not overpower the tree- level are interesting in their own right. However, 
in this section we return to the question of non-gaussianity in the cosmic microwave 
background. One can obtain a marginally tighter bound on by demanding that the 
loop expansion is stable for this expectation value as well as for the spectrum, and we 
shall see that increasingly stringent bounds are possible for higher correlation functions. 

The Schwinger formalism diagrams contributing to the three-point function are 
shown in Fig. [2j These diagrams make contributions to the three-point function which 
can be put into a form similar to Eqs. (14.71) and (I4.10j) for the two-point function. For 
example, for the (+, +, +) diagram we obtain 

( A^A^A^, 



(^ a (k 1 )5/(k 2 )^(k 3 ))* D 



x 



64 IL*? 

d 3 g d 3 r d 3 s 



A 3 



fi«JI(l + iA*77„)e" 



3^3 



P 3 (q, r, s)5(-ki - r + q)5(-k 2 - s + r)5(-k 3 - q + s) 



x Q + (k 1 ,k 2 ,k 3 ,q,r,s), (5.1) 
where now % G {1, 2, 3} and Q + is a function of the scalar momenta which is defined by 



fcl fc 2 &3 

ki+q + r k-2 - r + s fc 3 — q — s 



+ 



fcl fc)3 &2 

fcl + q + e fc 3 - q + s fc 2 - r - 



+ 
+ 



fc 2 fci k 3 

k'2 + r + s fci + q — r kz — q - 

k2 fc3 fci 

fc2 + r + s k:j + q — s fei — q - 



+ 
+ 



fc 3 fci fc 2 

k'j, + q + s fci — q + r ki — r — s 

fc3 fc2 fcl 

k:i + q + s fc 2 + r — s ki — q — T 



(5.2) 



The momentum polynomial P 3 (q, r,s) is the analogue of Eq. ( 14. 9 ft for the bispectrum, 
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Figure 2. Diagrams contributing to the three-point scalar expectation value. As 
before, straight lines (appearing on the external legs of the diagrams) indicate scalar 
quanta), whereas wavy lines indicate virtual gauge bosons circulating the the loop. 
These diagrams break into four groups of complex conjugate pairs, with the (+, +, +) 
and (—)—)—) diagrams forming one pair and the permutations of the (+,+,—) and 
(— , — , +) diagrams forming the other three groups. 



and depends symmetrically on each of the vector momenta q, r and s. Specifically, it 
obeys 

P 3 (q, r, s) = q 2 {r ■ s) 2 + r 2 (q ■ s) 2 + s 2 (q ■ r) 2 - (q ■ r)(q ■ s)(r ■ s), (5.3) 

whereas the time integrals can now be cast in the form of a six-parameter exponential 
integral, 



a± ct2 013 
Pi #2 P3 



c 



dr 



drj 



1 - ia lV )(l - ia 2 r) (1 - ia 3 C)e i/31 V^V^, (5.4) 



) C J— 00 1~ J— 00 v 

with the necessary deformations of the contour of integrations — or, more accurately, the 
inclusion of suppression factors in the exponentials — to cause convergence at large |^|, |t| 
or I ( I . We expect the same general considerations to govern the result of the integrals in 
Eq. ( 15. II) which applied for the two-point function: after integrating by parts, Eq. ( 15.41) 
can be expressed as an asymptotic series in inverse powers of 77*, together with nested 
integrals of the form J d( £ _1 e 1 ^ for some \x. 

We must also take account of the (+, +, )- and (— , +, +)-type diagrams. Since these 
are related by complex conjugation it suffices to consider only diagrams of one particular 
signature, for which we will choose the (+, +, — )-type, and then add in their complex 
conjugates. The — vertex can be attached to any of the external legs, corresponding to 
the any one of the momenta kj. The contribution from any such diagram, for example 
where the — vertex is associated with k 3 , can be written in a form similar to Eq. ( 15. ip . 
taking account of the fact that, irrespective of its numerical value, the time at the — 
vertex is taken to be later than the time of evaluation 7]*, and with Q + replaced by a 
different function Q_ which accounts for the necessary sign interchanges. We therefore 
arrive at 

(^(kO^M^Cks,)) d // " /A " A " V 



64 Ylik 
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x (1 + iA^)(l + ifc 2 t7„)(l - ifc 3 t7.)e~ M * x+ * a ~* !0 
d 3 gd 3 r d 3 s 



3^3 



-P 3 (p, q, r)5(-ki - r + q)<5(-k 2 - s + r)5(-k 3 - q + s) 



x Q_(k 1 ,k 2 ,k 3 ,q,r,s), (5.5) 

together with the equivalent expressions obtained by exchanging ki with k 2 and k 3 . In 
this expression, Q_ is defined by 

Q- = 



ki k2 —kz 

ki + q + r k-2 — r + s — kz — q — s 



+ 



&2 ki -kz 

k2 + t + s fci + q — r — &3 — g — s 



+ 
+ 



fel -k% k2 

ki + g + r — &3 — g — a k% — r + ; 

-kz fei fc 2 

— &3 — q — s k\ + q + r k2 — r + ; 



+ 



&2 — ^3 fel 

&2 + r + s —kz - q — s h>2 + q — r 

-kz k 2 fci 

— &3 — q — s k2 + r + s k\ + q — r 



(5.6) 



To see how this works in detail, it follows after integrating by parts that we can 
express Eq. (15.41) in the form 
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(5.7) 



where the coefficients £i to £9 are defined by 

£1 = i{9(a 2 - ai) + 3a 3 + 2(ft - 2ft) - ft} (5.8) 

^{9a 3 ft + 2ft 2 - 3a 3 ft - 4ft 2 - 2ftft - ft 2 + 3 ai (8a 3 + ft - ft - 5ft) 

(5.9) 



6 = 
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ft -ft 
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This function contains primitive fast divergences, but as in the case of the spectrum 
these will conspire to cancel among themselves in the final answer, leaving a result which 
contains only slow divergences which scale as powers of N*. As before, the combination 
fa+fa + fa is always independent of the loop momentum, whereas most of the integrals 
involve decaying exponentials of \fa\ or \fa\ and will play no role when the loop 
integral is dominated by its ultra-violet region. In Q + , there is no combination other 
than fa + fa + fa which can contribute, whereas in Q_ — for four of the terms — the 
combination fa + fa (and therefore fa on its own) also remains finite in the extreme 
ultra-violet limit. It follows that the terms involving £ 3 and can also contribute, 
although it will turn out that in Eq. (15.71) all these pieces cancel out of the final answer. 

In order to simplify the calculation from this point onwards, it is convenient to 
specialize to the equilateral limit in which all the ki are equal, giving k\ = &2 = k%. In 
practice this does not entail much loss of generality, since the approximations we are 
obliged to make in evaluating the time integrals in Eqs. (14. 8 j) and (15. 4p mean that we 
cannot take the ki to exit the horizon with a separation of more than a few e-folds. 
(Exactly similar remarks apply to the standard calculation of the tree-level three-point 
function.) One may now proceed by analogy with the scalar two-point function, as 
described in §H The three-point function is written in the form 

(5r(W(k 2 )5^(k 3 ))^ =k = (2vr) 3 5(^k l )^ rf ^3 f^s-J ^> ( 5 - 17 ) 

and the time- and momentum-dependence is absorbed into I3. This integral is somewhat 
more complicated than its counterpart for the two-point function. It is defined by 

h = Ak [ -^-P 3 (q,r,s){4A;(g + r + s) + (g 2 + r 2 + s 2 )ln27 
J q^r^s 6 L 

+ (qr + qs + rs)(4Af* + 4 7e + 9 In 3 - 4) - 2k 2 }, (5.18) 

with the specific assignments r = q — ki and s = q + k 3 . In addition, iV* satisfies the 
usual definition N* = \n\kr)*\, which is unambiguous here because of our assumption 
that the momentum triangle formed by the kj is equilateral. 

To proceed, we pick k x to point along the z axis of a spherical polar coordinate 
system. (We could also pick k 3 , but the result is independent of our choice.) In this 
system of coordinates the only non-trivial inner product we require is that of q with k 3 , 
which takes the form 

q ■ k 3 = fcgj cos^cos^is + sin 9 sin 813 cos(0 — 0i3)|, (5.19) 

where {0,4>} are the zenith and azimuth of q relative to ki, and {#i 3 ,0i 3 } are the 
corresponding zenith and azimuth of k 13 . Note that 0i 3 is devoid of significance, being 
only a coordinate choice, and cannot appear in any physical quantity, whereas 9 and 
cf) are simply variables of integration. The integral is quadratically divergent. With a 
cut-off A, we find 

h = -^A 2 (n. - \ hi 3) - l^A + % In -A> (4 + 2 cosfl 13 + V\ , (5.20) 
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where iV* is a slightly modified count of the number of e- folds since horizon exit, defined 



One might worry that the appearance of #13 in Eq. (15.201) represents of failure of 
symmetry between k 1; k 2 and k 3 . However, it should be remembered that once 
we specify any two of the momenta the third is determined by the triangle law 
hi + k 2 + k 3 = 0. Only the relative orientations of these vectors can enter physical 
quantities, and the absolute orientation of one vector and the azimuth of a second are 
simply coordinate choices. For an equilateral configuration, this leaves only one cosine 
which can play a role in physical quantities. 

The tree-level, single-field equilateral bispectrum is known to take the form [TTl [20] 



We drop power law divergences and neglect a possible threshold correction. This 
implies that, if we specialize to single-field inflation and return to the parametrization 
A(0) = exp(0/M ?i ), the one-loop corrected bispectrum can be written as 



where A' is again the proper scale of new physics, which can have a hierarchy with respect 
to the Planck scale of no more than around h^A'/if*) ~ 10. The power spectrum of the 
curvature perturbation, V-r, is computed using the tree-level scalar spectrum, which we 
assume is not destabilized by the loop correction from electric quanta. It follows that 
if, likewise, we do not wish the loop correction to be larger than the tree-level, we must 
have M" 1 < 50. 

This bound is somewhat stronger than is required to guarantee stability of the 
spectrum. Indeed, it is clear that successively higher-order expectation values will 
contain increasingly large negative powers of M^, whereas the leading correction to 
an n-point scalar expectation value comes from diagrams such as those in Fig. [2] where 
the electromagnetic quanta circulate in a ring between the external scalar legs. Such 
a diagram is always suppressed with respect to the tree-level by only one power of the 
loop-counting parameter (H/Mp) 2 , and therefore for very large n we will require 
to be closer and closer to unity. Does this imply that only = 1 makes sense for an 
effective field theory? In fact, this conclusion would be too strong. For too far below 
the Planck scale the present analysis does not show that the two- and three-point scalar 
expectation values are inconsistent with observation, but rather only that perturbation 
theory is insufficient to compute them. 

It is known that the tree-level power spectrum is an excellent match for observation, 
and that the temperature anisotropy is gaussian to high accuracy. It follows that 
there is a reasonably secure foundation to demand that the two-point function can 
be computed perturbatively. The situation for the bispectrum is less clear. It is known 



by 




(5.21) 




(5.22) 




(5.23) 
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that the non-linearity of the temperature anisotropy is not too large, although it may 
perhaps be of order / NL ~ 60 in the squeezed limit. This limit is the opposite of the 
configuration studied here, where one of the momenta is going to zero while the other 
two become approximately equal and opposite, but even in the equilateral case where 
estimation is more difficult there is a relatively stringent limit, |/nl| 300. There is 
not yet any hint from the data that the three-point function requires anything beyond 
conventional perturbation theory. Even less is known concerning the properties of the 
n-point expectation values for n larger than three. For this reason, it seems most 
conservative to adopt the approximate bound M^ 1 < 120 which corresponds to the 
applicability of perturbation theory for the power spectrum. 

6. Discussion 

In this paper I have computed the corrections to the spectrum and bispectrum of a 
set of light scalar fields which arise from coupling to an electromagnetic field during 
inflation. The calculation of inflationary observables is conventionally split into a 
quantum-mechanical initial condition, which specifies the momentum-dependence of 
expectation values at the time of horizon crossing, and a subsequent classical evolution. 
The interaction studied in this paper constitutes a correction to the initial condition, and 
can be thought of as arising from interactions with virtual quanta of the electromagnetic 
field. In a scenario where the gauge coupling is rolling rapidly during inflation, these 
virtual electromagnetic quanta experience strong amplification as they are drawn across 
the horizon and "freeze in" to the correlations among the scalar modes at that time. 

To which scenarios would this correction apply? The most direct application is to 
theories of so-called magnetogenesis, where one aims to produce a late-time magnetic 
field by breaking the conformal invariance of the U(l) gauge field. In this scenario the 
scalar fields which determine the expectation value A(</>) must be rolling on cosmological 
timescales, meaning that these scalars will typically also contribute to the curvature 
perturbation, (. It is £ which is visible as the adiabatic temperature perturbation in the 
Cosmic Microwave Background. Therefore, in order that the conventional perturbative 
predictions for the power spectrum and bispectrum are not destabilized we must demand 
that the correction from electromagnetic interactions does not dominate the tree-level. 
As shown in S@J this requires that M^ 1 < 120 in order that the spectrum is not 
destabilized, and the stronger condition (£J5]) M^ 1 < 50 for the bispectrum. Both of 
these conditions are approximately satisfied for conventional pictures of magnetogenesis, 
which, as shown in §|3l typically requires M^ 1 « 20. According to the discussion in §5] 
it is probably most conservative to adopt the limit M^ 1 < 120, because the evidence 
that current observation matches perturbative calculations of the bispectrum is rather 
less strong than for the spectrum. 

There is another mechanism by which synthesis of an electric or magnetic field could 
destabilize conventional perturbation theory. If the energy density in electromagnetic 
radiation begins to compete with the potential energy associated with the scalar 
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fields, then this would deform the comoving hypersurface on which inflation ends. 
This deformation would imply that £ was not conserved, but rather would acquire a 
contribution from this final hypersurface |73j. On the other hand, if the energy density 
in the electromagnetic field remains small, then we can be confident that there is no 
back-reaction on the metric. In this case, it follows that the hypersurface on which 
inflation ends suffers a negligible deformation and ( is given by the usual procedure, 
because the electromagnetic energy density never contributes to the expansion rate or 
the integrated number of e-foldings, N. 

When does energy in the electromagnetic field remain small? This question was 
addressed by Bamba [10] and Martin & Yokoyama [32], who gave a detailed discussion 
of the back-reaction problem. They concluded that when a magnetic field is synthesized 
there is never a problem with back reaction, whereas if the outcome is an electric field 
then the scenario is only free of instabilities associated with back reaction if the scale of 
inflation is taken to be very low. One way to understand this phenomenon is to note that 
the magnetic energy density ps ~ {dui) 2 is proportional to spatial gradients, whereas the 
electric energy density p# ~ cj 2 involves a time derivative. We must therefore include 
Pe at leading order in a gradient expansion of the scalar Klein-Gordon equation. If 
Pe becomes too large it begins to act as a source and causes an unwanted growing 
instability in the scalar fluctuations. In practice this "instability" may entail nothing 
more dramatic than a macroscopic flow of charge which shorts out the electric field [40J, 
but our ability to calculate is impaired and it is difficult to draw conclusions with any 
confidence. Therefore, unless the scale of inflation is taken to be very low, it seems that 
one should restrict attention to scenarios where A(0) grows during inflation. 

Couplings between scalars and U(l) gauge fields are invoked in certain theories 
of dark energy, such as the so-called "chameleon" [74] model. Such couplings were 
proposed in Ref. [51], following earlier work [75j [76], and studied using astrophysical 
means in Refs. [52], [53l [77]. (For laboratory limits, see Refs. [78], US]-) However, 
the chameleon is essentially inert during inflation: it rolls rapidly to its minimum, 
and is fixed there for the duration of inflation [80] . Since it is not in motion while 
observable scales are leaving the horizon, the parameter Q which controls the magnitude 
of the correction tends to zero. Moreover, because the energy density carried by the 
chameleon is cosmologically negligible, its expectation values do not contribute to those 
of the curvature perturbation. Therefore, although perturbation theory may fail for 
expectation values of the chameleon field if the associated mass scale is too small, this 
has no observational consequences. We conclude that no limit on the chameleon coupling 
(3 analogous to those obtained in Refs. (ED, [521 E3] can be extracted from the present 
analysis. 

Although fast instabilities were encountered at intermediate points in the 
calculation, the final expectation values were found to be free of fast divergences. 
Their time dependence was instead described by powers of N„, which measures by 
how many e-folds the fluctuation in question has passed outside the horizon. One can 
understand this as a consequence of a theorem due to Weinberg, which guarantees that 
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Figure Al. All-scalar loop correction from the V" vertex, common to all scalar field 
models of inflation. 



interactions involving gauge fields are well behaved in this sense |15j . In the present 
case, although the electromagnetic part of the interaction would be "dangerous" (in 
Weinberg's sense) as part of a scalar interaction because it involves time derivatives, it 
does not produce fast divergences. Indeed, the only effect of differentiation with respect 
to time when applied to a mode of the gauge field is to introduce factors of the slowly 
varying parameter Q. Such factors merely set the scale of the interaction, rather than 
changing the character of the time dependence, as would be the case for a scalar mode. 
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Appendix A. Scalar loop corrections from the V" vertex 

In this Appendix, I compute the scalar loop correction coming from a simple V" vertex. 
This loop correction is common to any inflationary model, and arises simply from the 
self-interactions implied by the scalar potential, unless V" is somehow tuned to be zero. 
This loop correction corresponds to the diagram in Fig. IA11 with two external scalar 
legs which are connected by a loop of circulating virtual scalar quanta. Just like the 
processes considered in the main text, this diagram comes in four different types labelled 
by the + or — flavours at each vertex, giving complex conjugate pairs (+,+), (— , — ) 
and (+,-), (-,+). 

The loop process described by Fig. lAll has already appeared in the literature, and 
has been the subject of detailed study by several previous authors, van der Meulen & 
Smit [62] studied this loop for both large and small internal momenta, finding corrections 
to the classical time dependence. These corrections are one possible source of the 
"quantum logarithms" described in the Introduction (31]). Bartolo et al. [25] considered 
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the same loop process, but focused on the infra-red region and found results in this 
regime which were in agreement with those of van der Meulen & Smit. In this Appendix 
the same loop is studied using the methodology applied in the main text, principally 
for the purpose of aiding comparison of the methods of the present paper with those of 
either van der Meulen & Smit or Bartolo et al. 

We work with a single-field model of inflation. The generalization to multiple fields 
is easy to obtain by identical methods, and keeping track of the necessary indices in 
the space of scalar fields only introduces needless notational clutter. It follows that the 
vertex is simple, and takes the form 



d3fci ( y fc3 W6(j2 k^50(k 1 )^(k 2 )^(k 3 ). 



(A.1) 



As with the two-point loop correction studied in §U it is sufficient to compute the 
contribution from the (+, +)- and (+, — )-type diagrams and then add in their complex 
conjugates. From the (+, +) diagram one obtains 

1 -,,A 2 f d 3 qd 3 r 



(50(k 1 )^(k 2 )}, D 



(K 



q3 r 3 



-<J(-ki + q-r)5(-k 2 



q 



X 



X 



;i + iA; 1 ^)(l + ifc 2 ^)e- ir '*^ +fc2 ) 

J fa — q — T kiq + far — qr fagr 
\ fa + q + T —k 2 q — far — qr k 2 qr ' 

whereas from the (+, — ) diagram one obtains 

d 3 g d 3 r 



fa — q — r fag + far — qr k 2 qr 
fa -f q + r —k±q — far — qr fagr 



(A.2) 
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-k 2 qr 
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fa + g + r —kiq — far — gr 
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(A.3) 



As in §|U these contributions are symmetric between q and r, and are both written in 
terms of a time integral which takes the form 



Qi o-i az 
Pi f3 2 As 



dr. 



iair] + a 2 r/ 2 + ia 3 77 3 )e iair? 



x(l 



iPiT + for 2 + i/3 3 r s )e^ T . (A.4) 

In order to evaluate the loop process, we shall need an expression for this integral in 
closed form. This can be obtained by iteratively integrating by parts, which produces a 
series containing "fast instabilities" associated with negative powers of rj, together with 
logarithmic contributions from integrals of the form f^* d( (~ 1 e ll3< '. This structure is the 
same as the one which we encountered when studying time integrals in §^fflSl although 
the result is slightly more complicated. This extra complexity stems from the fact that 
Eq. (1A.4|) contains more negative powers of r and rj than (for example) Eq. (14. lip and 
therefore requires more integrations by parts to separate its asymptotic dependence on 
77*, and also because the integrands contain more terms anyway. Fewer terms are present 
when dealing with the gauge field because of its conformal invariance, which gives the 



Magnetogenesis and the primordial non-gaussianity 



33 



associated gauge field wavefunctions a simpler temporal dependence than the minimally 
coupled scalar field. 

The time dependence which arises from Eq. (I A. 4 1) is somewhat complicated to 
express. One obtains 



«1 


a 2 


013 




[Pi 


$2 




"( 



o~i .o- 2 a 3 .a 4 cr 5 .cr 6 

s v* ??* ??* m V* 
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d.77 



(a lP3 + a 3 A 4 ) / — e iai " / — e 
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" dr , 
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(A.5) 



The parameters appearing in this expression involve certain combinations of the on and 
(3j, which are described by the choices 
1 



Ai 



\ 2 = p 1 X 1 ; A 3 = p 2 + AA 2 ; and A 4 = /3 3 + /3iA 2 



7i = Ai; 72 = A 2 + aiAi; 73 = A 3 + a 2 X\ - ctiA 2 ; 



74 = aiA 3 + a 2 X 2 — a 3 Ai; 75 = a 3 A 2 + « 2 A 3 

We also define a combination S, 

5 = at + Pi. 

In terms of these combinations, the pi satisfy 

cti(A 4 - pi) 



and 7 6 = a 3 A 2 



A 4 

Pi = y; 



P2 



P3 = "2A 4 + «ip 2 . 



We now have all the ingredients necessary to define the af. 



01 



7i, 
6 ' 



0-9 



72 - 5(Ti 
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a 4 = 



7 4 - 8a 3 + pi 



0-5 



73 + 5a 2 
4 



75 + 8°~4 ~ P2 



and 



o"6 = 7e + Sa 5 - p 3 . 



(A.6) 
(A.7) 
(A.8) 

(A.9) 
(A.10) 

(A.ll) 

(A.12) 
(A.13) 



Let us focus first on the infra-red region, where the momentum circulating in the 
loop is small compared to the momenta ki = k 2 = |k| on the external legs. In this limit, 
we can take q ~ and r « fc. Although Eq. (1A.5I) contains a large number of primitive 
fast divergences (up to and including in the pure power-law sector, and up to and 
including r]~ 2 when multiplied by additional logarithmic divergences, bearing in mind 
that purely imaginary divergences cancel out of the final answer), these all cancel in 
both the ultra-violet and infra-red. The leading time dependence instead comes from 
the double exponential integral. This can be written 



J -00 V J 



_ el /3ir = ln 2| ai ^| + 7E + ln 

T Z 



a i 



Pi 



In 1^*1 +/(«!/&), (A.14) 



where f(x) is a complicated function of its argument which can be given in closed form 
in terms of the Euler-Mascheroni constant when x — 1. We will focus only on the 
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leading divergence as 77* — > 0. Keeping only this term, one finds that the loop corrected 
power spectrum takes the form 

f 2 fV"'\ 2 N 2 \ 
p ( k ) = PoM h + g [^) jf 2 ^k£j, (A.15) 

where Po(k) is the tree-level power spectrum and I is an infra-red regulator which has 
been introduced to prevent a divergence in the momentum integral. As described in 
Refs. (EH E21 E31 EH EH [25] , this regulator can be understood as the size of a comoving 
"box" in which we perform the quantum field theory calculation. This box must be 
chosen to be sufficiently large that the region of space for which we wish to obtain a 
prediction can fit comfortable inside it, but small enough that all fluctuations remain 
perturbatively small. Eq. ( 1A.15I) is identical with Eq. (1) of Ref. [25J. 

In the opposite limit, where very hard virtual quanta are circulating in the loop, q 
and r are both much larger than k. In this limit, one finds that the loop correction has 
the form 

P,^,^^™^. (A, 6 ) 
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